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Abstract 



We reformulate the dual gauge model of the long-distance Yang- 
Mills theory in terms of the two-point Wightman functions formal- 
^ . ism. In the flux-tube scheme of Abelian dominance and monopole 

condensation, the analytic expressions of both monopole- and dual 
gauge boson-fields propagators are obtained. Formally, new features 
of higher order quadratic equations for a monopole field are given. 
Finally, we show how the rising-type potentials in the static system of 
color charges are naturally derived. 

PACS 12.38.Aw, ll.15.Kc, 12.38Aw, 12.38Lg, 12.39Mk, 12.39Pn 

1 Introduction 

There is a general statement that the color confinement is supported by 
the idea that the vacuum of quantum Yang-Mills (Y-M) theory is realized 
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by a condensate of monopole-antimonople pairs [1]. In such a vacuum the 
interacting field between two colored sources located in x\ and x<i is squeezed 
into a tube whose energy E tube ~ \x\ — x 2 \. This is a complete dual analogy 
to the magnetic monopole confinement in the type II superconductor. Since 
there is no monopoles as classical solutions with finite energy in a pure Y-M 
theory it has been suggested by 't Hooft [2] to go into the Abelian projection 
where the gauge group SU(2) is broken by a suitable gauge condition to its 
(may be maximal) Abelian subgroup U(l). It is proposed that the interplay 
between a quark- ant iquark p£iir is cLiicilcL gtous to the interaction between a 
monopole-antimonopole pair in a superconductor. 

In fact, the topology of Y-M SU(N) manifold and that of its Abelian 
subgroup [[/(l)]^" 1 are different, and since any such gauge is singular, one 
might introduce the string by performing the singular gauge transformation 
with an Abelian gauge field [3] 

A^x) -»• A^x) + -f d^{x) , (1.1) 

where Q(x) is the solid angle subtended by the closed space-like curve de- 
scribed by the string at any point x = (x°, x 1 ), and g = 2n/e is responsible for 
the magnetic flux inside the string, e being the Y-M coupling constant. Here, 
we choose a single string in the two-dimensional (2d) world sheet y x (r, a), for 
simplicity. Obviously, the Abelian field-strength tensor F^ v = d^A u — d u A fl 
transforms as 

F*(x)^F*(x)+G, u (x) , 

where a new term 

G yv (x) = -$-[d lt ,d v \n(x) , 

is valid on the world sheet only [4] 

G^(x) = I e^ uaP J J dadr ' V ^ 5 4 [x - y(a, r)] . 

Formally, a gauge group element, which transforms a generic SU(N) connec- 
tion onto the gauge fixing surface in the space of connections, is not regular 
everywhere in spacetime. The projected (or transformed) connections con- 
tain topological singularities (or defects). Such a singular transformation 
( |1 . 1| ) may form the worldline(s) of magnetic monopoles. Hence, this sin- 
gularity leads to the monopole current J™ on . This is a natural way of the 
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transformation from the Y-M theory to a model dealing with Abelian fields. 
A dual string is nothing but a formal idealization of a magnetic flux tube 
in the equilibrium against the pressure of surrounding superfluid (Higgs-like 
field) which it displaces [5,6]. 

Recent lattice results [7] give the promised picture that the monopole 
degrees of freedom can indeed form a condensate responsible for the confine- 
ment. The expression for the static heavy quark potential, using an effective 
dual Ginzburg-Landau model [8], has been presented in [9]. In the paper 
[10], an analytic approximation to the dual field propagator without sources 
and in the presence of quark sources, and an expression for the static quark- 
antiquark potential were established. 

The aim of this paper is to consider the model in 4d based on the dual 
description of a long-distance Y-M (LDY-M) theory which shows some kind 
of confinement. We study the model of Lagrangian where the fundamental 
variables are an octet of dual potentials coupled minimally to three octets 
of monopole (Higgs) fields. The dual gauge model is studied at the lowest 
order of the perturbative series using the canonical quantization. The basic 
manifestation of the model is that it generates the equations of motion where 
one of them for the scalar Higgs field looks like as a dipole-like field equa- 
tion. The monopole fields obeying such an equation are classified by their 
two-point Wightman functions (TPWF). In the classical level there is some 
intersection with the Froissart model [11] containing the scalar field satisfy- 
ing the equation of the fourth order. In the scheme presented in this work, 
the flux distribution in the tubes formed between two heavy color charges 
is understood via the following statement: the Abelian monopoles are ex- 
cluded from the string region while the Abelian electric flux is squeezed into 
the string region. 

In Sec. 2, we introduce the essence of the dual gauge Higgs model and 
the classical solutions. In our model there are the dual gauge field C£(x) and 
the monopole field (i = 1, ...,N C (N C — l)/2; a=l,...,8 is a color index 

) which are relevant modes for infrared behaviour. The local coupling of the 
Bi-field to the C^-field provides the mass of the dual field and, hence, a dual 
Meissner effect. Although C^{x) is invariant under the local transformation 
of [/(l)^- 1 C SU(N C ), = ■ H is an SC/(Ay -gauge dependent object 
and does not appear in the real world alone (N c is the number of colors and 

— * 

H stands for the Cartan superalgebra) . The commutation relations, TPWF 
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and Green's functions as well-defined distributions in the space 5(3? ) of 
complex Schwartz test functions on §l d will be defined in Sec. 3. In Sec. 
4, we study the monopole- and dual gauge- field propagations. In Sec. 5, 
we obtain the asymptotic transverse behaviour of both the dual gauge field 
and the color-electric field. The analytic expression for the static potential 
is obtained in Sec. 6. Sec. 7 contains the discussion and conclusions. 



2 A dual Higgs gauge model and classical so- 
lutions 

The dual description of the LDY-M theory is simply understood by switching 
on the dual gauge field C^(x) and the three scalar octets Bi(x) (necessary to 
give mass to all and carrying color magnetic charge) in the Lagrangian 
density (LD) L [12] 



L = 2Tr 

where 



W BA , (2.1} 



F^v = d^C u - d u C^ - ig[C^ C v \ , 
D ll B i = d li B i -ig[C li ,B i ] , 

and Bi are the SU(3) matrices, g is the gauge coupling constant of the 
dual theory. The Higgs fields develop their vacuum expectation values (v.e.v.) 
Bq. and the Higgs potential W{Bj) has a minimum at B 0i . The v.e.v. B . 
produce a color monopole generating current confining the electric color flux. 
It is known, the LD ( |2.1| ) can generate classical equations of motion carrying 
a unit of the Z3 flux confined in a narrow tube along the z-axis (corresponding 
to quark sources at z = ±00). This is a dual analogy to the Abrikosov [13] 
magnetic vortex solution. 

As the next step we introduce the color ansatz 

^ = Ec;k, (2.2) 

a Z 

where the vector potential is dual to an ordinary vector potential in the 
Y-M theory, (l/2)A a are generators of SU(3). Following paper [12] in the 
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sense of representing the quark sources by the Dirac string tensor G^ iv (x) 
having the same color structure as in (|2.2|) , one can arrive at a more suitable 
form of the LD Q) 

L{G, V ) = ~G%, + 4|(3„ - igC,)<j>\ 2 + 2(<9 M 3 ) 2 - W{<j>, 3 ) , (2.3) 
where 

G^u = d^C v — d u C^ + G^ u , 

and 4>(x) and 03 (x) denote the complex scalar monopole fields. We choose the 
color structure for the QCD-monopole field Bi (belonging to the fundamental 
representation of SU C (3)) like in [12], and the effective potential stands 

W(B,B,B 3 ) = ^A{11[2(5 2 + B 2 - Bl) 2 + (B 2 3 - B 2 ) 2 ] 

+1\2(B 2 + B 2 )+Bl-W 2 } 2 } , 

where 

= 01 = 02 = B 1>2 - lBl,2 , 03 = B 3 , 

while A provides the weak couplings between the scalar fields. 
The dual gauge field C M satisfies the relation 

d,C u - d u C, =* (d,A u - 8 V A,) 

in the absence of charges, and the duality transformation is realized by in- 
terchanging the gluon field A^(x) and C^(x). In the maximally Abelian 
gauge [14] An(x) = A a ^(x){T a /2) behaves as the Abelian gauge field C M (x) = 
A^(x)(t 3 /2) approximately because the off-diagonal fields are suppressed by 
the gauge transformation. The LD ( [2.3| ) is invariant under the local gauge 
transformation of the dual gauge field 

C^x) -> C^x) - d^A(x) 

and the phase transformation of the QCD-monopole field 

01,20*0 exp[-z gA(x)} 0i |2 O*O , 

where A(x) is the real field in S(^t 3 ) at any fixed x°. The local gauge sym- 
metry is spontaneously broken because of (0) o ^ in ( |2.3| ). The generating 
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current of ( [2.3| ) is nothing but the monopole current confining the electric 
color flux 



J^T = Aig[<f>\d ll -igC^(t>- (Kdp + igCJP 
which enters into the equation of motion in the form 

3 



d u G^{x) 



J: 



x 



jmon 



0, means that 



The formal consequence of the J™ on conservation, d^J r u 
monopole currents form closed loops. 

To find a solution of this model, one can consider the monopole-field as 
a solution where (B(x)) Q ^ 0, (B(x)) ^ 0, (B 3 (x)) ^ 0. We choose 

B(x) = b(x) + B ,B(x) = b(x) ,B 3 (x) = b 3 (x) + B 

with the boundary conditions at large distances p from the center of the flux 
tube with $ as an angle in the cylindrical coordinates [12] 



C — > e# 

99 



B exp (ifi) ,B 3 ^B as p — > oo . 



(2.4) 



In terms of the fields b(x), b(x), b 3 (x) and C^x) the LD fl2.3p is divided into 
two parts 



Li + L 2 



(2.5) 



where L\ in the lowest order of g and A and with the minimal weak interaction 
looks like 



-C 2 
3 ^ 



1 



W + {d.bf + -(d,b 3 ) 



+ m 2 C 2 - -p 2 (50b 2 + 181Q + Smd^bC^ . 



(2.6) 



Here, m = gB and p = v2\Bo are masses of the dual gauge field and the 
monopole field, respectively. The remaining part of (|2.5| ) turns out to be 



L 2 = 8g (d^b)C,b - (d^b)C,b + Ag 2 (d^C u - d v Crf ■ (b 2 + b 2 + 2B b) 
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-^A[25(6 4 + b A ) + 96* + 100B & • (b 2 + b 2 + B 2 ) + 28B (b 2 + b 2 + bb 3 ) ■ b 3 

+36B b 3 (bl + Bl) + 26 2 (256 2 + 7bj) - 2^(506 + 186 3 ) + 56B 2 bb 3 + Ub 2 b 2 3 ]. 

Let us consider the canonical quantization of ( j2.6|) . The equations of motion 
are 

(A 2 + /i 2 ) b(x) = ; 



A 2 b(x) +m(d-C) = ; (2.7) 

(A 2 + f i 2 2 )b 3 (x) = 0; 

(A 2 + m\) C,(x) - d,(d ■ C) + 12mdJ>(x) - d v G^{x) = , (2.8) 

where fi 2 = (50/3) /z 2 , /i 2 = 12 /i 2 ,m 2 = 3 m 2 . By taking the divergence of 
(|2.8|) one can get 

A 2 b{x) - ^-d^G^x) = , (2.9) 
9m 

while the formal solution of equation Q2.8| ) looks like 

C ll (x) = a&'G ia ,{x)-pd ll b{x) , 

where a = (3 m 2 ) -1 , (5 = 4/m. We see that the dual gauge field is defined via 
the divergence of the scalar field b(x) shifted by the divergence of the Dirac 
string tensor G fJ/l/ (x). At the same time, we propose in the standard manner 
the Dirac string which can be understood as a straight line connecting two 
objects with opposite charges. For large enough x, approaching such a string, 
the monopole field is going to its v.e.v. while C^(x — > oo) — > 0. Hence, 

J™ on {x^oc) ^8m 2 C^ . 
For a very weak C^-field one can say that in the d — 2h dimensions 

A 2h b(x)~0, h = 2,3,..., (2.10) 
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but 

A 2 b{x) ^0 . 



Here, the solutions of equation ( 2.10|) obey locality, Poincare covariance and 



spectral conditions, and look like the dipole "ghosts" at h=2. In this model, 
the role of the dipole field at d=4 is held in the d=2 dimensions by the 
simple pole field, and the analogy of the behaviour between the d=2 and 
d=4 dimensions can be found at the level of Wightman functions in the free 
case, at least [15]. This is related to the model proposed in [16] in a study of 
the Higgs phenomenon, from which the present model is distinguished by the 
coupling to the dual gauge field C^x), and the gauge field strengh tensor is 
shifted by the Dirac string tensor G^(x). Thus, we obtain that the massless 
scalar field b(x) occurs in the model since the symmetry realizes such a way 
that the LD (|2.3| ) is invariant under the local gauge transformations above 
mentioned but (b(x)) Q ^ 0. 

Our aim is to find Green's function of a scalar field b(x) obeying Eq. 
( |2.10| ). The propagator r h (x) is defined via TPWF Wh(x) in the d=2h- 



dimensions 

r h (x) = (Tb(x) 6(0) > = 6(x°) W h (x) + 6(-x°) W h (-x) , (2.11) 

where 

W h (x) = (b(x)b(0)) (2.12) 

is the distribution in the Schwartz space S'(3l 2h ) of temperate distributions 
on ^R, 2h and obeys the equation 

A 2h W h {x) = . (2.13) 

The general solution of eq. ( |2.13|) should be Lorentz invariant and is given 
in the form [17,18,16] at h=2 

W 2 (x) = a x In / + a \ + a 3 , (2.14) 
—x,, + % e x u x, — % e x u 

where (i=l,2,3) are the coefficients to be defined later, while / is an arbi- 
trary parameter with dimension minus one in mass units. Its origin becomes 
more transparent from [15, 19]. In fact, the first and the second terms in 



8 



( |2.14| ) are related to the scalar dipole field and the scalar pole field, respec- 
tively [15]. In the first case, the solution is 



Wo(x) 



-i E 2 (x) 



(4tt) 



In 



—x 2 + i e x° 

H 1 



while for the second term the solution looks like 

1 1 



WoX) 



-iD 2 (x) 



(2 



7r Y —xf, + i ex 







(2.15) 



(2.16) 



3 TPWF as classical distributions 

Before going into the quantization procedure, let us briefly consider the clas- 
sical distributions (in the sense of generalized functions [20]) of TPWF and 
T-ordered TPWF for the 6(x)-field at large x\. The TPWF ( gT2|) at h=2 is 
provided by the distribution 8(p°) 5'(p 2 ) as [21] 

W 2 (x) ~ J d 4 p 6(p°)5'(p 2 ) exp(-tpx) ~ - ln[-/i 2 x 2 + i6(x )] 

= - ~\n\fj, 2 x 2 \ + insgn(x°)6(x 2 )} , (3.1) 

where jx is an arbitrary parameter (the infrared regularization parameter). 
The distribution 9{p°) S'(p 2 ) in ( |3.1|) is uniquely defined only on the test 
functions f(p) G S o (&0 (S (^) = {/(p) G S(^), f(p = 0) = 0}) 



J d 4 p 6(p°) 5\p 2 ) f(p) = hm A / d 4 p 9{p°) 5\p 2 - v 2 ) f(p) , 
d 3 p 1 



r+ 2p° 2 (np) 



(np) 



(nd) 



and for an arbitrary fixed time-like unit vector n„ in the case we choose 



n,, 



[1,0) from V + where 



x G $1 : x > \x\ < £ 



1/2 
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is an open upper light cone in the M-space. Under the dilatation transfor- 
mation x — > ax (a > 0) the TPWF (|3.1|) acquires the additional term 



Wo(x) 



W 2 (ax) = W 2 {x) 



1 



2 (2 tt) 



In a. 



(3.2) 



It could be interpreted as a spontaneous symmetry breaking of the dilatation 
invariance of ( |2.10| ). This is an important point in the special role of the field 
b(x). 

In general, the r^o^-function fl2.11|) is a well-defined distribution in S'^R, 2h ) 
and obeys the /i-ordered quadratic differential equation in the d-dimension 
of space-time 

A 2h r h d (x) = S d (x) , h = 2,3,... , 



A 



Q2 Q2 
+ ... + 



(J 2 



dx 2 i ' ' dx 2 

in the following cases [20,22] 
i) even d and h > d/2 



dx 2 ,. 



m+l 



dx 2 



m+n 



Th (X) 



vrf/2-1 



exp[(7r/2) in] 



A h (h - d/2)! (h - 1)! 



Jx P(x) + i e 



-d/2+h 



■In 



fi P(x) + ie 



(3.3) 



where P(x / 
ii) even d and h < d/2 



,-yi 2 I I ,-v> 2 ™ 2 

O/ ^ \ • • • \ <m ^ 



r h d (x) 



h exp[(?r/2)in] 



/} P(x) + z e 



rl 2" ft , 



(3.4) 



In the case of d = 2/i-dimensions, expression ( |3.3|) can be more simplified as 
follows: 



v ; v ; A h (h-1)\ 



li 2 P(x) +ie 



For the following nearly realistic consideration of the confinement occurrence, 
we are interested in the case i) ( [3. 3D which allows one to study highly sin- 
gular objects for the confinement-like picture. In view of that, for h — 2 in 
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the M-space following Zwanziger [23] the Fourier inversion F[t 2 m (x)} of the 
distribution t 2 m (x) looks like 



T2 M (X) 

and obeying the equation 



7T 2 8 



2 dp 1 



p M ln(/ 2 p 2 



(p 2 + % e) 



(3.5) 



T2 M (X) 



1, P 2 = ijpf-^ =1 p 3 ^ 



where / = exp(7e — 1/2) • (2/i) 1 (je stands for the Euler constant). Note 
that in terms of weak derivatives the distribution (|3.5|) can be rewritten as 
[15] 



ThHx) 



lim 



(_l)d/2 



(p 2 - k 2 + i e) d l 2 



+ in d/2 5 d {p) In 



K 



4/i 



/~/2 



for any h allowed in the dimension d. The remaining case ii) 
the following singular object: 



(3.6) 
leads to 



r h d {x) 



lim u~ d (-X 2 + 



I K 



where /x 2 X 2 = p\ + ... + p 2 m — p 2 m+1 — ... — p^ +n . In the space M one has 



Th M (x) 



lim u 



r.-4 



A 2 



— p 2 + ie 



4 A dual field propagator 



For calculation of the coefficients a\ and a 2 in Q2.14Q one has to introduce 
invariant functions E 2 (x) = E 2 (x) — E 2 (—x) and D 2 (x) = D 2 (x) — D 2 {— x) 
[23] (see also formulae (|2.15| ) and (|2.16|) ) which define the commutator written 
in the general form [16] 



b{x), 6(0) = (2 tt) 2 i [4 ai E 2 {x) + a 2 D 2 {x)} 



(4.1) 



E 2 (x) = (8n)- L sgn(x u )6(x 2 ) , 
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D 2 (x) = (2 7T)- 1 sgn(x°)S(x 2 ) . 
In the space S"(!R 4 ) on 3ft 4 the propagator for the 6(x)-field looks like [16] 



T 2(x) = a,\ In | /x 2 x^ | + iir9(x 



a 2 



x 2 + iir5(x 2 ) 



a 3 . 



(4.2) 



The coefficients a\ and a 2 in Q4.2|) can be fixed using the canonical commu- 
tation relations (CCR) 

^(0)] | = i^5 3 (f) , 



6(X),7T 5 (0) 



i<5 3 (£) 



respectively. Here, the conjugate momenta 7Tc Ax) and vr^(x) look like 



n- b (x) = 8 +mC°(x) 



4 

The direct calculation leads to (see also Appendix): 

,2 



(4.3) 



m 



a i 



a 2 



12(2tt) 2 
1 

~6(2tt) 2 



Restricting the W 2 (x) function ( |2.14| ) to only the first term, one can obtain 
that r 2 (x) satisfies the equation 



A 2 ) 2 fAx) = i5Ax) 



(4.4) 



where f 2 (x) = 3(2/m) 2 t 2 (x). The formal Fourier transformation in S"(3ft4) 
gives 



f 2 (p) = weak lim 

k 2 <<i 6 (2 ir)^ 



m 



1 



K 



(p 2 — k 2 + i e) : 



+ i7r In — 5Ajp) 
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Here, k is a parameter of representation and not the analogue of the infrared 

we used the well-known mathematical 



mass, k 2 = K 2 /p 2 . To derive 
trick with the generalized functions [20,16,15] 

weak lim ! a^p expHp* 

ft 2 «l J 



:-v, 



(p 2 — k 2 + it) 



2 i 



T(h) (47T) 



-{In 2 - lE - \h{kJ-xI + i e) 



+0[-k x^-k x^ ln(K yj-x* +ie)]} . 

To go into the structure of the dual C^-field propagator, we need to define 
the general form of the commutation relation [C M (x), C„(y) ]. To do this 
procedure let us consider the canonical conjugate pair {C^, nc v } (f° r the LD 
where 7Tc m (x) results from ( f4.3| ). The consequent CCR looks like 



z^5 3 (f),(4.6) 



-C^x), d v C (0) - d o C u (0) - g v(d ■ C(0)) + A 0u {0) 



where A 



X 



9p,v (d- C(x)) — G^vix) tends to zero as x — > and the Dirac 



string tensor G fll/ (x) obeys the equation 

(A 2 + M 2 ) G^(x) = , as x 



Here, the equations of motion ( f2.7| ) and ( |2.9D were used, and M = 3m. 
Obviously, the following form of the free C^-field commutator (see Appendix) 



; C M (x), 0,(0)] = ig^ £ ml E 2 (x) + c D 2 (x) 



(4.7) 



ensures the CCR 



at large x 2 . Here, both £ and c are real arbitrary 



numbers. The TPWF for C^x) stands as (see also [15,19]) 

w^ u {x) = (C^(x)C u (0)) Q = ig^ v ^m\E 2 {x) + cD 2 (x) + c 1 F 2 (x) + c 2 ,(4.8) 

where F 2 (x) = ( a(x) a(0) ) = ix 2 is TPWF for the harmonic field a(x), c\ 
and c 2 are arbitrary real numbers. In fact, the third and the fourth terms 
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in ( |4.8| ) did not appear in ( [4.7| ) due to the properties of the F 2 ~ (^-function, 
namely, F 2 (x) = F 2 ~(x) — F 2 {—x) = and the triviality, respectively. One 
can easily verify that ( |4.6[ ) and ( |4.7j ) can be used to derive the following 
requirement for £ 



4 



(4.9) 



The free dual gauge field propagator in S"(9ft 4 ) in any local covariant gauge 
can be assumed as 



= / d 4 x exp(ipx)r^ u (x) 



C 7 P 2 + i e 



^ ml kip) + ct 2 (p) 



(4.10) 



where 



\je ) 



(4tt) 2 



£ m\ ln(-/i 2 x 2 + z e) + 



xl + ie 



tiijp) = weak lim 

K 2 <<1 



1 



(p 2 — k 2 + i e) 2 



+ i 7r 2 In 



Hp) 



t 2 (p) = weak lim - 



s 2 «i 2 (p 2 - k 2 + z e) 

The gauge parameter ( in ( |4. 10|) is a real number. According to Eq. 
the following requirement on Green's function t^x) 



(A 2 ) 2 f fiu (x)=t5 4 (x) , 

leads to that a constant c has to be equal to zero where f^ u (x) = r^ u (x) / (£ m 2 ). 
At the end of this section, let us consider the v.e.v. 



V»(x) = ([j™ n (x),b(0)}) 



(4.11) 



where the monopole current expressed in terms of the 6-field and G^ v Dirac 
tensor is 



- d^A 2 b(x) -a A 2 VG^x) - d»G^{x) 



m 



14 



Using the commutation relation ( |4.1| ) and the well-known relation between 



the following distributions like A 2 sgn(x°) 



9(x 2 ) = 4sgn(x°) S(x 2 



we arrive 



at the formal expression for V^(x) 



%m _ 
o u 
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sgn(x°) S(x 2 ] 



(a A 2 - V 



9"G^(x),6(0) \ .(4.12) 



It is easily seen that the requirement of the Goldstone theorem for occurrence 
of a S^pfJ term in the Fourier transformation of is satisfied due to the 
presence of the first term in ( |4.12|) which gives the Fourier transformed term 
~ p^sgn(p°) 5(p 2 ). The second term in ft4.12|) comes from the dual string 
effect. 



5 The LDY-M solution for a gauge field 

In this section, we consider the approximate topological solution for the dual 
gauge field in the LDY-M theory. The equations of motion are 

d u G, u = 6zg[<f ) *(d,-tgC,)<P - <P(d, + z g C,) <f>*] ; (5.1) 

(fy - i 9 C,) 2 4> = ~ A (32 B 2 - 25 | | 2 - 7 2 ) . (5.2) 

At large distances Eq. (|5.2|) transforms into the following one: 

50 

( d ,-zgC,) 2 <f ) = Y X(B 2 - |0| 2 )0. 

To get the solution, let us do the polar decomposition of the monopole field 
0(x) using new scalar variables x( x ) an d f{%) 

4>(x) = -j= exp(z/(x)) [x(x) + B ] . 

Then, the equation of motion ( |5.1|) transforms into the following one: 

d»G^ = 6g( X + B ) 2 (gC,-d,f) . (5.3) 

This means that the 6(x)-field is nothing but a mathematical realization of 
the "massive" phase B ■ f(x) at large enough distancese: 

b( x ) = ^S(x)f(x)+ 
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m f [2S{x)-l]C li (x)dx li + -^— f [A 2 C^x) - d"d^C u (x)]dx^ 1 , (5.4) 
where S(x) = (1 + x{ x ) / Bq) 2 . Now we have to define the flux as 



n = y G^ u (x)da^ , (5.5) 

where a^ v is the 2d surface element in the M-space. A substitution of C^(x) 
from ( |5.3D into (|5.5| ) gives 



n= f i 5_1(a:) + / G^ u (x)da^ + -^ d^f(x)dx» ,(5.6) 



where T means the large closed loop where the current d^C v — d v C^ is can- 
celed. Integrating out over the loop V in the third term in ( |5.6| ) is nothing 
but the requirement that the phase f(x) is varied by 2 7cn for any integer 
number n associated with the topological charge [24] inside the flux tube. 
One can present the fields C^x) and <f>{x) in the cylindrical symmetry case 
using the radial coordinate r (see also [24]) 

C — > eg , — > (r) , 



and f = 2n n = n6 with 9 being the azimuth around the z-axis. Thus, the 
field equation looks like 



d 2 C(r) 1 dC(r) 



3m 2 [3 + 2S(r)]C(r) + QnmB S(r) = . (5.7) 



dr 2 r dr 
The following boundary conditions: 

4 n i— 
C(r) = — , X (r) = (V2-l)B , asr^oo; 

C(r) = , x( r ) = , as r — > 

are conjugate with those in (p.4[) . At large enough r >> /i" 1 ((// = v^2A"i?o) _1 
defines the transverse dimension(s) of a monopole field around the tube) Eq. 
transforms into the following one: 

d 2 C(r) 1 dC(r) , . „ d2 



dr 2 r <ir 



+ 3<?(4n-7<7 (7)^ = 
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with the asymptotic transverse behaviour of its solution 

C(r) ~ J — — e kmr 1 + — ; , fc = V21 . 

The field C (r) grows rapidly when the radial distance from the center of the 
flux tube r < r$ ~ 3 fm and approaches An/ 7 g as soon as r > tq. Now 
it will be very instructive to clarify the singular properties of the field b(x) 
related to the phase f(x) by means of ( |5.4j ). Since the phase is provided by 
the asimuth 6 around the z-axis, V/ = (n/r) e$. In the strong limit of the 
mass m of the C M field (as well as the mass of the monopole field) at large 
distances 

V x Vb = 2irnB 5(x)5(y)e z , 

where e z is the unit vector along the axis z and the 5-functions stand for 
the center of the flux tube [24]. Hence, the real singular character of the 
field b(x) is confirmed by its singular behaviour at the center of the flux tube 
for nonzero monopole condensate. At the end of this section we give the 
transverse distribution in r of the color electric field E z {r) (see also [24]) 

E c = VxC=- e z = E z (r) ■ e z 

r ar 

which has the following profile in the flux tube at large r 




6 Static potential 

In this section, we intend to obtain the confinement potential in an analytic 
form for the system of interacting static test charges of quark and antiquark. 
Our statement is based on the dual character of the field C^(x) to a gluon 
field where C^(x) is just the interacting field provided by the monopole field 
b(x) and the divergence of G^ u (x). We have found that b(x) plays the role 
of the dipole-type field at h — 2 (see Eq. (|2.10|) ). The mass of the dual field 
Cfj,(x) is nonzero and equal to m = g B . It is assumed that the mass m 
results from the Higgs-like mechanism when the dual field C M (x) interacts 
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with the <t>(x) field, namely an octet of dual potentials C M coupled weakly 
with three octets of scalar fields (f) a (x)(a = 1, 2, 3). 

According to the distribution ( |3.6j ), the first term in the expansion for 
the static potential 



J d 3 pe^ F{r h \x)} Wo 



in 3? 3 is a rising function with r = \x\ [25,22] 

Pstatir) ~ T(3/2 - h) r 2 ^ 3 . (6.1) 

It is obvious that at h > 2 the distribution ( |6.1| ) increases with r linearly 
(h = 2) or faster {h > 2). Let us represent the distribution r°" as the Taylor 
series around some regularization point ctq 



1 + cj ln(/xr) + -cu 2 ln 2 (/ir) + ... , (6.2) 
2 



where uj = a — a is an infinitesimal positive interval at a ^ —d, —d — 2, 
The potential (|6.1j) is simplified to [22] 



P s tat(r) ~ ^-r[l + ln(/ir) + ...] . 

O 7T 

The Fourier transform F{r a } of ( |6.2|) into S'(dtd) for the limit u — > leads 
to the following singular distribution in the whole region of the existence of 
the analytic function r a at o ^ —d, —d — 2, ...: 

11 \P 2 ) r(-a/2) ' 

In general, the static potential is defined as 

P atat {r)= lim L(r) , (6.3) 
where the action A(r) is given by the colour source- current part of LD 

L{p) = -±{-p)T, v {p)t{p) ■ 



18 



It is known that for such a system of heavy particles the sources are given 
by a c-number current 

j£(x) = Q a g" [5 3 (x - xi) - 5 3 (x - x 2 )\ 

with Q a = e p a being the Abelian color-electric charge of a quark while 
p a is the weight vector of the SU(3) algebra: p\ = (1/2, v3/6), p 2 = 
(-1/2, \/3/6), p 3 = (0, -1/V3) [24]. Here, £ a and a?2 are the position vectors 
of a quark and an antiquark, respectively; the label a=(l,2,3) corresponds to 
the color electric charge. The calculation of the potential ( |6.3|) is most easily 
performed by taking into account the Fourier transformed quark current 



J Ma (p) = 2 tt Q Q *(pP) (e-** - e^ 2 ) , (6.4) 
and by making use of the representation in the sense of generalized functions 



weak lim 

K 2 <<1 



(p 2 — k 2 + i e) 



1 d 2 



A dp 2 —p 2 — it 



In 



-p —is 



2 (p 2 + i e) 



+ in 2 In — 5 4 (p) 
-p 2 - k' 



5 - 3 In 



. (6.5) 



Due to the presence of the <5(p°)-function (see (|6.4| ) in evaluating of the action 
A(r), the remaining 3-dimensional integral over d 3 p will be easily calculated 
out using the instructive formula [20,15] 



d 3 pe ipx p m In? p 



r(/? + i)(-i) 1 



d 



„ m +3 £ T (P - i + l)r(i + 1) V^m, 



H m Wr, (6.6) 



where 



2 m+3 ^3/2 



r[(m + 3)/2] 



m ^ —3, —5, ... 



T(-m/2) 

r = \x\, p = |p|, x G 3? 3 , p € 9ft 3 . 

As a consequence of the dual field propagator ( (4.1Q ) the static potential 
(|6.3|) at large distances looks like 



Pstat(r) 



16 71 



£m 2 r[5 + 6(A + lnpr)] + (^\\ 



(6.7) 
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where A = 1 + -0F- (1/2 + 2 ^7^) - (2^+ 1) In 2 < and the last term 



in ( |6.7|) is just the positive correction and not the analogue of the Coulomb 
part of the potential due to the one-gluon exchange. Neglecting the last term 
in ( |6.7|) and taking into account formula (|4.9| ), one can conclude 

Pstat(r) ~ |^m 2 r(-12.4 + 6 ln/ir) . (6.8) 

64 7T 

Hence, the string tension a in P s t a t{^) = ar emerges as 

a ~ — m 2 (-12.4 + 3 In At ) , jl>9m, (6.9) 



16 \ m 



where r in the logarithmic function in ( |6.8|) has been changed by the char- 
acteristic length r c ~ 1/m which determines the transverse dimension of the 
dual field concentration, while jl is associated with the " coherent length" in- 
verse and the dual field mass m defines the "penetration depth" in the type 
II superconductor. For a typical value of the electroweak scale [i ~ 250 GeV 
we get a ~ 0.20 GeV 2 for the mass of the dual C^-field m = 0.6 GeV and 
a s = e 2 /(47r)=0.37 obtained from fitting the heavy quark-antiquark pair 
spectrum [26]. The value of the string tension ( |6.9| ) is quite close to a phe- 
nomeno logical one (eg., coming from the Regge slope of the hadrons). Making 
the formal comparison of the result obtained here in the analytic form let us 
remind the analogue with the well-known expression of the energy per unit 
length of the vortex in the type II superconductor [27,9] 

(i = t^ 1 J^Y (6.10) 

32 n 2 \m A J 

where 0o is the magnetic flux of the vortex, and are penetration 
depth and the coherent length inverse, respectively. On the other hand, the 
string tension in Nambu's paper (see the first ref. in [1]) is given by 

8 7r y m v z J 

with m s and m v being the masses of scalar and vector fields and g m is a 
magnetic-type charge. It is clear from formula (|6.7|) that for a sufficiently 
long string r >> mT x the ~ r-behaviour of the static potential is dominant; 
for a short string r << mT x the singular interaction provided by the second 
term in ( |6.7D becomes important if the average size of the monopole is even 
smaller. 
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7 Summary and discussion 



We studied the dual gauge model of the long-distance Yang-Mills theory in 
terms of two-point Wightman functions. The quantization of the model has 
been provided by using CCR, thus avoiding other methods (e.g., path in- 
tegral use). We intended to give our understanding of the confinement by 
making use of nothing else but the well-known tools of quantum field theory 
based on LD given in [12] as well as the renormalization model and symme- 
try properties. Among the physicists dealing with the models of interplay 
of a scalar (monopole, Higgs) field with a dual vector (gauge) boson field, 
where the vacuum state of the quantum Y-M theory is realized by a con- 
densate of the monopole-antimonopole pairs, there is a strong belief that the 
flux-tube solution explains the scenarios of color confinement. Based on the 
flux-tube scheme approach of Abelian dominance and monopole condensa- 
tion, we have obtained the analytic expressions for both the monopole and 
dual gauge boson field propagators ( |4.5| ) and ( |4.10| ), respectively, in S"(9? 4 ). 
These propagators lead to a consistent perturbative expansion of Green's 
functions. However, the Fourier transformation of the first term in TPWF 
( |2.14|) gives the occurrence of the 5' (^-function. This is a consequence of 
the nonunitarity of the translations, and the spectral function with such a 
term gives an indefinite metric [16]. We have found that the Goldstone the- 
orem is valid in our model in the form taking into account the Dirac's string 
effect ( |4.12j ). In fact, we obtained that the characteristic 5{pfy term natu- 
rally appears in the Fourier transformation of the v.e.v. of the commutator 
( 4.11 ) of the monopole current and the scalar local field b. In principle, a 
similar result has to be expected if one replaces the 6-field in (|4.11| ) by any 
product of the gauge field C M and b. We see that the fields b(x) and b 3 (x) 
receive their masses and the b(x) field in combination with d u G llv {x) form 
the vector field C^(x) obeying the equation of motion for the massive vector 
field with the mass m = g B . The solution of the 6(x)-field can be identified 
clS db "ghost" -like particle in the substitute manner. 

The occurrence of an arbitrary parameter I in (|2.14|) and ( |2.15|) leads 
to breaking their covariance under the dilatation transformation ( |3.2|) and 
provides spontaneous symmetry breaking of the dilatation invariance of Eq. 



( 2.10 ). The monopole condensation, formulated in the framework of LDY-M 
model, causes the strong and long-range interplay between heavy quark and 
antiquark, which gives the confining force, through the dual Higgs mecha- 
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nism. We have obtained the analytic expression for the static potential (|6.7|) 
at large distances. The form of this potential grows linearly with the distance 
r apart from logarithmic correction. The latter comes from the second term 
in the expression ( |6.5|) (see also flSHD ). 

Making an analytic comparison of ei fl6.10|) and e 2 ( |6.11| ) with a in 



one can conclude that we have obtained a similar behaviour of the string 
tension a to those in the magnetic flux picture of the vortex and in the 
Nambu scheme, respectively, as well as in the dual Ginzburg-Landau model 
[9]. 

Finally, it is to be noted that we have played the game with the choice 
of the gauge group where the Abelian group appears as a subgroup of the 
full Y-M gauge group. This is a very instructive method of calculating the 
confinement potential in the static limit in the analytic form. However, we 
understand that no real physics can depend on such a choice. Now, there is 
a next step in more formal consideration of the Y-M theory where it seems 
to be a new mechanism of confinement [28,29]. 
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9 Appendix 

The following relations [21,16,15] were used for calculation of the coefficients 
ai and a 2 in (14. II): 



D 2 (x° = 0, x) = , dpD 2 (x). =g ^S 3 (x); 
(A 2 ) 2 E 2 (x) = d 2 E 2 (x) ] = d E 2 (x) ] = E 2 (0,x) = 

la; u — la; u — 

c%E 2 {x\ o = 8ng 0fl 5 3 {x) , A 2 E 2 (x) = D 2 {x) . 
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